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In  order  to  make  future  satellites  both  smaller  and  smarter,  more  navigation  information  must  be  extracted 
from  simpler,  smaller  sensors.  One  of  the  simplest  sensors  is  an  optical  or  infrared  camera.  With  a  camera,  a 
satellite  can  track  a  second  satellite  located  within  its  field-of-view.  This  simple  measurement  is  the  foundation 
of  angles-only  navigation.  By  its  very  nature,  angles  only  navigation  cannot  determine  the  relative  range  to  an 
object.  Even  as  the  dynamics  associated  with  orbital  rendezvous  and  proximity  operations  unfold,  the  relative 
range  will  generally  remain  unobservable.  In  this  paper  we  confirm  that  an  angles  only  navigation  system 
can  observe  range  if  small  maneuvers  can  be  executed,  and  we  show  that  the  level  of  accelerometer  accuracy 
determines  how  well  the  range  can  be  observed. 


Schmidt  1  of  20  23rd  Annual  AIAA/USU 

Conference  on  Small  Satellites 
American  Institute  of  Aeronautics  and  Astronautics  Paper  SSC09-VI-3 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

AUG  2009 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2009  to  00-00-2009 

4.  TITLE  AND  SUBTITLE 

Improving  Angles-Only  Navigation  Performance  by  Selecting  Sufficiently 
Accurate  Accelerometers 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Utah  State  University, 1400  Old  Main  Hill, Logan, UT ,84322 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

SSC09- VI-3,  AIAA/USU  Conference  on  Small  Satellites,  Logan,  Utah,  August  2009 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

- 1 - 

- 1 - 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 
ABSTRACT 

Same  as 
Report  (SAR) 


18.  NUMBER 
OF  PAGES 

20 


19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


I.  Nomenclature 

Accents,  subscripts  and  superscripts 

True  Value  =  b 
Estimated  Value  =  3 
Measured  Value  =  o 
Previous  Time  step  Estimate  =  o~ 

Updated  Estimate  =  6>+ 

Chaser  camera  frame  =  ocam 
Chaser  body  frame  =  oh 
Chaser  values  =  oc 
Target  values  =  ot 
Cross  Product  Matrix  Form  =[ox] 

Errors  and  noise 

Measurement  Error  (v)  =  o  —  d 
Residual  Error  (e)  =  o  —  o 
Misalignment  (e)  =  vector  of  small  angles 
Bias  ( (3 )  =  vector 

Process  Noise  (w)  =  continuous  white  noise 
Measurement  Noise  (77)  =  continuous  white  noise 

II.  Introduction 

The  purpose  of  this  study  is  to  show  how  accelerom¬ 
eter  accuracy  affects  the  performance  of  a  Kalman  fil¬ 
ter  processing  angles-only  measurements  during  orbital 
rendezvous.  While  angle  measurements  can  be  obtained 
from  other  types  of  sensors  (Lidar  or  Radar  for  example), 
the  camera  has  the  additional  advantage  of  being  entirely 
passive. 

By  its  very  nature,  an  angle  measurement  contains  no 
range  information.  If  angle  measurements  are  made  to 
multiple  objects  with  known  locations  that  lie  in  differ¬ 
ent  directions  this  limitation  can  be  overcome.  When  the 


angle  measurements  are  made  to  a  single  object,  such  as 
during  orbital  rendezvous,  the  lack  of  range  information 
can  become  a  serious  limitation. 

Various  strategies  have  been  developed  to  overcome  this 
limitation.  Two  of  the  most  successful  techniques  are 
taking  “apparent  diameter”  measurements,  and  perform¬ 
ing  translational  maneuvers.1  Both  methods  have  their 
disadvantages.  Apparent  diameter  measurements  are 
range  limited  based  on  the  resolution  of  the  camera  and 
require  prior  knowledge  of  the  satellite  being  observed. 
Translational  maneuvers  consume  fuel  and  lower  the 
life  of  the  satellite.  In  this  paper  the  translational  burn 
method  will  be  considered. 

Lingering  Questions  on  Angles -Only  Navigation 

Angles-only  navigation  (AON)  has  been  researched  for 
use  on  spacecraft  since  the  days  of  Gemini,2  and  proved 
useful  for  orbital  rendezvous  during  the  Apollo  years.  It 
was  successfully  implemented  autonomously  on  Deep 
Space  I.3  AON  is  useful  because  is  does  not  require 
heavy  and  expensive  sensors  like  radar.  Instead,  an  ex¬ 
pensive  radar  is  replaced  by  a  video  camera  and  occa¬ 
sional  maneuvers. 

However,  there  are  still  serious  questions  about  AON  use 
during  orbital  rendezvous.  For  distances  greater  than 
100km,  Gauss’  method  may  be  an  option  to  determine 
range.4  However,  at  closer  distances,  Gauss’  method 
breaks  down.  This  is  because  whole  families  of  trajecto¬ 
ries  will  exhibit  nearly  identical  line-of-sight  (LOS)  mea¬ 
surement  histories  as  seen  in  figure  1 .  They  only  differ 
in  their  range  component.  If  this  motion  is  linearized  us¬ 
ing  CW  equations,  then  the  LOS  measurement  histories 
will  truly  be  identical.  Without  a  unique  solution  for  the 
given  measurements,  the  relative  state  of  the  chaser  will 
remain  unobservable. 

The  behavior  of  Kalman  filters  when  processing  LOS 
measurements  has  been  analyzed  well  studied.1,5  It  has 
been  shown  that  the  component  of  the  covariance  ma¬ 
trix  parallel  to  the  LOS  vector  will  grow  almost  with¬ 
out  bound  unless  maneuvers  that  change  future  measure¬ 
ments  are  performed.  Within  a  short  time,  state  estima¬ 
tion  can  return  bogus  results  (See  figure  8). 


Schmidt  2  of  20  23rd  Annual  AIAA/USU 

Conference  on  Small  Satellites 
American  Institute  of  Aeronautics  and  Astronautics  Paper  SSC09-VI-3 


Figure  1.  A  family  of  relative-motion  trajectories  exhibiting  identical  line-of-sight  measurement  histories. 


Observability  Burns 

In  order  to  estimate  range,  satellites  using  AON  for  or¬ 
bital  rendezvous  must  make  maneuvers  in  order  to  ob¬ 
serve  relative  range.  As  Figure  2  shows,  only  one  range 
will  satisfy  the  known  change  in  position  <5r  resulting 
from  a  known  acceleration  delivered  by  the  spacecraft’s 
thrusters  as  long  as  dr  is  not  parallel  to  rnominai.  The 
change  in  the  LOS  vector  for  a  given  Sr  is  known  as  the 
observability  angle  (0)6  .  These  LOS  measurements  are 
processed  in  an  EKF  for  state  estimation. 

The  difference  between  the  actual  6  and  the  calculated  6 
is  a  function  of  accelerometer  and  image  centroiding  er¬ 
rors.  As  shown  in  Figure  3,  if  the  calculated  observability 
angle  is  too  small,  then  the  Kalman  filter  may  tend  to  ig¬ 
nore  it,  because  it  falls  within  noise  of  the  system.  Thus, 
accelerometer  performance  may  have  a  tremendous  im¬ 
pact  on  the  ability  of  the  navigation  filter  to  estimate  rel¬ 
ative  range. 

Kalman  Filter 

The  Kalman  filter  was  first  developed  by  R.E.  Kalman7 
and  has  been  improved  and  expanded  over  the  years  to 


process  all  sorts  of  measurements.  A  Kalman  filter  es¬ 
sentially  propagates  an  estimated  state  (x)  and  the  covari¬ 
ance  ( Px )  of  that  state  in  real  time.  When  measurements 
are  made,  the  estimated  state  is  improved  and  the  size  of 
the  covariance  decreases. 

When  angle  measurements  are  processed  in  a  Kalman 
filter,  no  information  can  be  gleaned  parallel  to  the  LOS 
vector.  Applying  the  filter  in  a  rendezvous  situation  with 
no  observability  burns  will  result  in  a  very  long,  skinny 
covariance  ellipse.  As  a  result,  extremely  large  and  very 
small  values  will  be  contained  in  the  same  covariance 
matrix.  Within  a  short  time  numerical  issues  can  crop  up 
during  the  measurement  update  portion  of  the  Kalman 
filter  algorithm. 

To  help  overcome  this  issue,  the  Kalman  filter  will  be 
reformulated  to  update  the  square  root  of  the  covari¬ 
ance  matrix.  This  formulation  is  known  as  a  square-root 
Kalman  filter,  and  has  the  significant  advantage  of  mak¬ 
ing  large  numbers  smaller  and  small  numbers  larger — 
making  the  filter  much  more  robust  when  propagating  a 
long,  skinny  covariance  ellipse.8  The  specific  formula¬ 
tion  the  square  root  extended  Kalman  filter  (SREKF)  is 
adapted  from  Tapley9  (See  figure  4). 
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Figure  2.  Change  in  observability  angle  as  a  result  of  dr  produced  by  a  known  acceleration 
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Figure  3.  Observability  angle  as  a  function  of  time 


III.  Methods 

A  satellite  rendezvous  simulation  created  in  Mat- 
lab/Simulink  is  used.  It  includes  noise  and  bias  on  sen¬ 
sors  and  actuators.  Vehicle  dynamics,  environmental 
models,  sensor  models,  and  actuator  models  already  ex¬ 
ists  as  Simulink  blocks.  Only  the  navigation  algorithms 
was  developed  for  this  research. 

Essentials  include  a  satellite  with  three  sen- 
sors-accelerometers,  line-of- sight  camera,  and  star- 
camera.  Actuators  include  thrusters  and  momentum 
wheels.  This  satellite  is  modeled  with  familiar  Euler 
equation  rigid-body  dynamics  and  J2  gravity  with  noise. 
The  flight  software  guidance  is  waypoint  driven  for  trans¬ 


lation  and  target  tracking  for  attitude.  Attitude  control  is 
obtained  with  a  phase-plane  controller  for  thrusters  and 
a  PID  controller  for  momentum  Wheels.  Translational 
control  leverages  a  PD  controller  for  station  keeping,  and 
Clohessy- Wiltshire  (CW)  equations  targeting  for  trans¬ 
fers.  Navigation  is  detailed  in  Section  III.A. 

III.A.  Kalman  Filter  Development 

A  Square-Root  EKF  is  used  for  navigation.  The  oper¬ 
ations  of  the  square-root  filter  is  distinguished  from  the 
standard  filter  as  seen  in  figure  4.  The  filter  has  been  for¬ 
mulated  to  process  line-of-sight  measurements.  The  fil¬ 
ter  assumes  J2  gravitational  dynamics  and  uses  onboard 
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accelerometers  and  star-tracker  measurements  to  directly  locity  of  both  the  chaser  and  a  target  vehicle  that  it  is 
propagate  the  states.  The  filter  runs  onboard  a  chaser  tracking.  The  filter  will  also  estimate  bias  and  misalign- 
spacecraft  and  will  estimate  the  inertial  position  and  ve-  ment  of  the  accelerometer  and  cameras. 


Figure  4.  Square-Root  Extended  Kalman  Filter  (SREKF)  and  Extended  Kalman  Filter  (EKF)  flow  charts  (See  section  III.B) 
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III. A.  1.  Filter  Design  Model 

This  is  the  navigation  system  designer’s  ’’working 
model.”  This  model  accounts  for  the  position  and  veloc¬ 
ity  of  the  chaser  and  target  vehicles,  misalignment  and 
measurement  noise  on  the  accelerometers  and  the  optical 
camera,  and  bias  on  the  accelerometers.  It  uses  process 
noise  in  the  vehicle  acceleration  channel  to  account  for 
unmodeled  effects  like  drag  and  solar  pressure. 

The  only  measurements  that  will  be  processed  in  the 


Kalman  filter  directly  will  be  LOS  and  accelerometer 
measurements.  Accelerometers  will  be  used  to  propagate 
position  and  velocity  states  directly.  As  a  result,  the  ac¬ 
celerometer  measurement  noise  (rjc)  will  be  treated  like 
a  process  noise. 

The  filter  design  model  can  be  written  as: 


x  =  f(x)  +Bw 

y  =  h(%)  +  1 1  cam 


rc 

Vc 

vc 
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(2) 


The  measurement  equation  is  detailed  below.  See  figure 
5  for  an  explanation  of  the  azimuth  and  elevation  angle 

values.  Ti^b{q)  =  (5) 


h(x)  = 


1 _ 

j^cam  jj^cam 

tan(el) 

j^cam  jj^cam 

RreP  ~  1/3x3  [^camx]](Tl^b(rt  rc)  rcam)  ~ 
\j^cam  j^cam  j^cam^T 


+  <?1  -  <?2  “  ^  2(4142 -40*73)  2(^143+4042) 

2(4142  +  4043)  4o  —  4i  +  42  —  43  2(4243-4041) 

(3)  [  2(4143-4042)  2(4243  +  4041)  4o~  <71  “42+43 

(4) 


where 


INERTIAL-TO-BODY  TRANSFORMATION  MATRIX 
The  transformation  matrix  comes  directly  from  the  star 
camera,  which  returns  a  quaternion  measurement  ( q)\ 


qo 

qi 

qi 


L  <73 


4x1 


cos(6/2) 
usin(  |) 


u  =  unit  vector  defining  axis  of  rotation 


0  =  angle  of  rotation  in  radians 


(6) 
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Figure  5.  Azimuth  and  Elevation  Measurements  in  Chaser  Body  Frame 


State  and  State  Covariance  Propagation 
Equations  These  are  the  equations  the  filter  will  use 
to  propagate  the  estimated  states  and  covariance 


k  =  f{x) 

/>+!  =  gjggT  +BQB'8t 


(7) 


The  strength  of  the  process  noise  Q  is  related  to  the  pro¬ 
cess  noise  w  by  the  relationship: 


Q8{t'  -t)=E[w{t')w(t)T] 


(8) 


Where  the  E  operator  is  the  expected  value,  and  8(t'  —  t ) 
is  the  Dirac  delta  function. 
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Q  may  be  represented  as  a  matrix  as  shown  below. 


Q  = 


(9) 

The  values  in  these  submatrices  usually  come  from  the 
hardware  specification,  while  Q ^  and  Q ^  are  func¬ 
tions  of  unmodeled  accelerations  of  the  chaser  and  tar¬ 
get  which  would  would  include  J?>+  gravity  effects,  drag, 
solar  radiation  forces,  etc. 

The  linearization  of  the  dynamics  equation  is  needed  to 
solve  for  the  transition  matrix  <£.  The  system  dynamics 
equation  (/(*))  is  linearized  as  follows: 
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Each  element  of  the  above  matrix  is  a  column  of  partials  as  seen  below. 
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Many  of  these  partial  derivatives  are  zero,  resulting  in  the  following. 
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The  individual  partial  derivatives  are  evaluated  as  follows: 
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In  the  Kalman  filter  the  estimated  states  (x)  will  be  used 
to  evaluate  these  partial  derivatives. 

Once  F  is  computed,  the  transition  matrix  is  given  by: 


<t>  =  eFdt  =  /  +  Fdt  - 


F2dt 2  F3dt 3 


2! 


3! 


(15) 


III.B.  State  and  State  Covariance  Update  Equation 

The  update  equation  is  where  the  standard  and  square- 
root  Kalman  filters  differ.  A  comparison  of  the  two  up¬ 
date  algorithms  may  be  seen  in  figure  4. 


The  standard  update  equation  works  well  on  long 
wordlength  machines,  but  can  cause  the  covariance  ma¬ 
trix  to  lose  positive  definite  nature  during  the  update 
step  when  the  covariance  matrix  is  ill-conditioned.8  The 
square-root  Kalman  filter  update  equations  are  consid¬ 
erably  more  complex  then  the  standard  update  equa¬ 
tions,  but  they  improve  numerical  accuracy  and  main¬ 
tain  the  positive  definite  nature  of  the  covariance  ma¬ 
trix.  The  square  root  algorithm  used  in  this  paper 
was  adopted  from  Tapley  and  Maybeck.8,9  The  ma¬ 
trix  square-root  was  calculated  by  way  of  Cholesky  de¬ 
composition,  which  is  referred  to  in  figure  4  as  “chol.” 
Cholesky  decomposition  results  in  a  lower  or  upper  tri- 
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angular  matrix  where  each  is  the  complex  transpose  of 
the  other.  For  the  lower  triangular  case,  the  result  times 
its  transpose  results  in  the  original  matrix.  Cholesky  de¬ 
composition  is  among  the  most  stable  of  matrix  opera¬ 
tions.10 


velocities.  So  the  measurement  sensitivity  matrix  be¬ 
comes: 


H  = 


dh(x) 


0  0  0 


dSi 

d  Beam 


2x21 

(19) 


The  individual  partials  are: 


R  =  chol(P) 

RRT  =P 


(16) 

(17) 


III.B.l.  Linearization  of  Measurement  Equation 

The  linearization  of  the  measurement  equation  results 
in  the  measurement  sensitivity  matrix  (H),  which  is  re¬ 
quired  to  solve  for  the  Kalman  gain  (K). 


dh(x) 

drc 

ajt, 

drt 

^rel 

u  cram 


Rl  ”  Rx 
_?z  _L  () 

-  Rx  R*  -I  2x3 

1/3x3  \fcam  X  ]]  7/—^  / 

1/3x3  — 

[(r^bir,  -  rc)  -  rbcaf  x] 


(20) 


Note  that  the  H  matrix  will  take  the  general  form: 


dh(x) 

~  d*  x~ 

dh{x)  dRbrd 

dRbrel  i 

(18) 

dh(x) 

dJk 

^k  ^k 

dRbrel 

d_$d 

dRbre, 

drc 

dvc 

drt  dvt 

dP 

deacc 

d^cam. 

The  angle  measurement  is  not  a  function  of  acceleration 
bias  (/?),  accelerometer  misalignment  (eflCC),  or  vehicle 


IV.  Model  Implementation 

The  simulation  tool  of  choice  is  Matlab/Simulink.  See 
Figure  6  for  an  example  of  this  simulation.  The  values  of 
he  variables  in  the  equations  in  section  III  are  found  in 
abies  1  through  5.  It  is  important  to  note  that  key  noise 
parameters  and  initial  covariances  must  be  matched  as 
good  as  possible  to  the  “real”  values  in  the  simulation  in 
order  for  the  filter  to  perform  as  expected. 


ExternalData 


Figure  6.  Example  of  simulation  model  implemented  in  Simulink. 


IV.A.  Noise  and  Time  Constants  chaser  was  varied  in  order  to  produce  the  results  found 

in  section  V. 

The  noise  strengths  and  time  constants  are  the  nominal 
values.  The  accelerometer  measurement  noise  on  the 
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Table  1.  Process  Noise,  Measurement  Noise,  and  ECRV  time  contents 


Parameter 

Description 

Value  in  Simulation 

Value  in  Filter 

Units 

Q(Oc 

Strength  of  Process  Noise  due  to 
random  accelerations  on  Chaser 

0 

10"18 

km 2 
s 3 

Qric 

Strength  of  Process  Noise  due  to 
accelerometer  measurement  noise 
being  processed  directly  in 
propagator 

10“12 

10“12 

km 2 

P 

QcOt 

Strength  of  Process  noise  due  to 
random  accelerations  on  Target 

0 

10-18 

km 2 
s3 

to 

O 

O 

"CD 

Strength  of  Process  noise  on 
accelerometer  bias 

10“20 

10-20 

km 2 

P 

Q  &)acc 

Strength  of  Process  noise  on 
accelerometer  misalignment 

0 

o 

<N 

1 

O 

rad 2 
s 

Q®cam 

Strength  of  Process  noise  on 
camera  misalignment 

0 

5  x  10-13 

rad2 

s 

R 

Variance  of  Measurement  Noise  on 
Angle  Measurements 

10“6 

10“6  0 

0  10“6 

rad 2 

^accfi 

ECRV  time  constant  for 
accelerometer  bias 

106 

106 

sec 

^ acc 

ECRV  time  constant  for 
accelerometer  misalignment 

oo 

106 

sec 

% cam 

ECRV  time  constant  for  camera 
misalignment 

oo 

106 

sec 

IV.B.  Initial  Conditions  ble  2  dictates  the  value  of  the  initial  covariance  of  that 

variable. 

Note  that  the  components  of  the  initial  conditions  are  all 
3x1  vectors.  The  value  of  the  standard  deviation  in  ta- 


Table  2.  State  Initial  Condition 


State 

Description 

Standard  Deviation 

Value  in  Filter 

Units 

fc 

Chaser  Position  in  ECI  Frame 

10-2 

II 

X  ^ 

5114.067 

-3998.013 

-335.012 

km 

Vc 

Chaser  Velocity  in  ECI  Frame 

10“5 

II 

X  ^ 

4.8199 

6.1714 

0.07119 

km/s 

ft 

Target  Position  in  ECI  Frame 

10-2 

II 

X  ^ 

5114.3258 

-3997.682 

-335.0156 

km 

Vt 

Target  Velocity  in  ECI  Frame 

10“5 

^  Vi  X 

II 

4.81954 

6.171718 
_  -.071167  _ 

km/s 

p 

Bias  on  accelerometers 

10-20 

0 

km/s1 

£acc 

Misalignment  on  accelerometers 

10-1° 

0 

rad 

^cam 

Misalignment  on  camera 

10“4 

0 

rad 
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Table  3.  Covariance  Initial  Conditions 


Table  4.  Navigation  filter  constants 


Component  of 
Covariance 

Value 

Units 

p 

rrcrc 

(10-2)2x/3,3 

km 2 

PvcVc 

(10-5)2x/3x3 

km2 /s2 

Pr,r, 

(10  ~2)2xI3x3 

km2 

PvtVf 

(10 -5)2x/&3 

km2 /s2 

Ppp 

10-20x/3x3 

km2 /s4 

Pc  p 

^acc  t-acc 

10-20x/3x3 

rad 2 

p 

P-cam  P-cam 

10-8  x  /3*3 

rad 2 

Symbol 

Definition 

Constant  and  Units 

Gravitation  Constant 

398600.4415  km’ /s1 

h 

Second  order 
gravitation  parameter 

0.0010826269 

Re 

Radius  of  the  Earth 

6378.1367  km 

n 

Constant  for 
calculations 

[0  0  1] 

dt 

Propagation  step  size 

0.25  sec 

Table  5.  Navigation  filter  inputs 


Symbol 

Definition 

Units 

a 

measured  acceleration  from 

accelerometers 

( km/s 2) 

q 

measured  attitude  from 

star  camera 

unit  quaternion 

y 

Angle  measurements  to 

Target  Satellite 

unitless  2D  vector 

7W 

transformation  matrix  (body  to  inertial) 

See  section 

V.  Analysis  and  Results 

V.A.  Observability  Maneuvers  Performance 

In  order  to  judge  the  effect  of  accelerometer  accuracy  on 
Kalman  filter  performance,  the  simulation  was  run  for 
the  flight  path  seen  in  figure  7.  In  brief,  the  satellite  per¬ 
forms  no  maneuvers  for  the  first  500  seconds.  After  this 
it  burns  positively  and  negatively  in  the  crosstrack  di¬ 
rection  successively  every  75  seconds.  The  importance 
of  these  maneuvers  is  emphasized  by  examining  figures 
8  and  9,  which  shows  navigation  performance  with  and 
without  accelerometers. 

The  effects  of  using  poor,  average,  and  good  accelerom¬ 
eters  on  navigation  performance  is  shown  if  figures  10 
through  12  and  may  be  contrasted  with  the  “perfect  ac¬ 
celerometer”  results  found  in  figure  9.  These  results  were 
generated  with  a  square-root  Kalman  filter  with  all  values 
defined  as  seen  in  section  IV,  except  for  the  accelerome¬ 
ter  accuracies,  which  are  defined  in  table  6. 

A  couple  of  other  methods  were  investigated  as  well. 
The  first  is  to  improve  state  estimation  by  increasing  the 


magnitude  of  the  thrust  maneuvers  so  that  they  are  more 
easily  distinguished  from  accelerometer  noise  and  error. 
The  second  method  is  to  ignore  the  accelerometer  mea¬ 
surements  when  a  burn  is  not  being  performed.  These 
two  cases  are  shown  in  figures  13  and  14  for  average  ac¬ 
celerometers. 

The  final  plot,  figure  15,  shows  results  when  LN-200 
type  accelerometers  are  used.  These  commercial,  off- 
the-shelf  accelerometers  have  a  acceleration  noise  vari¬ 
ance  of  3  x  10~13  km2 /s4,  slightly  better  than  the  av¬ 
erage  case,  and  a  bias  with  a  standard  deviation  of 
3  x  10 ~6  km/ s2,  much  worse  than  the  10 ~20  km /s2  used 
for  all  the  other  runs.11 


Variance  of  Acceleration 

Measurement 

Perfect 

0  km2  /s4 

Good 

10  ~14  km2 /s4 

Average 

10  ~l2km2/s4 

Poor 

10  -wkm2/sA 

Table  6.  Poor,  average,  and  good  accelerometers  specifications 
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Figure  7.  Flight  path  for  accelerometer  comparison 


3-o  Covariance  and  error  vs  time 
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Figure  8.  Relative  position  covariance  and  errors  with  no  maneuvers 
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;elerometer  Measurements  (km/s  )  body  z-axis  (vert)  (km)  body  y-axis  (horiz)(km) 


Figure  9.  Relative  position  covariance  and  errors  with  perfect  accelerometers 
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;elerometer  Measurements  (km/s  ) 


Figure  10.  Relative  position  covariance  and  errors  with  good  accelerometers 
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body  x-axis  (boresight)  (km) 


3-o  Covariance  and  error  vs  time 
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Figure  11.  Relative  position  covariance  and  errors  with  average  accelerometers 
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Figure  12.  Relative  position  covariance  and  errors  with  poor  accelerometers 
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;elerometer  Measurements  (km/s  )  body  z-axis  (vert)  (km)  body  y-axis  (horiz)(km) 
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Figure  13.  Relative  Position  Covariance  and  errors  with  higher  thrust  levels 
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;elerometer  Measurements  (km/s  ) 
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Figure  14.  Relative  Position  Covariance  and  errors  with  active  accelerometers  only  when  thrusting 
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3-o  Covariance  and  error  vs  time 
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Figure  15.  Relative  Position  Covariance  and  errors  with  LN-200  type  accelerometers 


VI.  Conclusions 

As  expected,  the  results  show  a  correlation  between  nav¬ 
igation  performance  and  accelerometer  accuracy.  If  the 
acceleration  level  of  the  maneuvers  are  not  overwhelmed 
by  the  measurement  noise  then  the  filter  can  estimate 
range.  If  the  acceleration  levels  are  overwhelmed  by  the 
the  noise,  then  the  filter  will  fail  to  estimate  the  relative 
state  accurately  (See  figure  12). 

It  is  common  to  improve  filter  performance  by  process¬ 
ing  the  accelerometer  measurements  only  when  a  burn 


is  being  commanded.  Increasing  the  magnitude  of  the 
thrust  is  also  effective  at  improving  navigation  perfor¬ 
mance.  Comparing  figures  11,13  and  14  reveals  that  pro¬ 
cessing  accelerometer  data  only  when  maneuvers  are  be¬ 
ing  executed  is  slightly  more  effective  at  improving  state 
estimation  then  increasing  burn  strength  and  both  tech¬ 
niques  are  better  then  the  original  case  (figure  11). 

The  run  with  LN-200  type  accelerometers  (figure  15) 
also  converges  to  the  correct  state  despite  the  increased 
bias  on  the  accelerometers. 

These  results  are  preliminary  and  further  testing  and 
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investigation  will  be  required  before  the  effects  of  ac¬ 
celerometer  error  and  processing  on  state  estimation  can 
be  fully  quantified.  However,  these  results  indicate  that  a 
camera  combined  with  average  commercial  off-the-shelf 
accelerometers  and  occasional  maneuvers  may  be  able  to 
estimate  the  relative  state  accurately  enough  for  proxim¬ 
ity  operations,  and  deserves  further  investigation. 
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